close all

%%
xx=0:dx:alx;
yy=0:dy:aly+dy;
load([saldir,'initial.mat'])
t=0;tx=['time=',num2str(t)];
figure('name',tx)
pcolor(yy,xx,squeeze(wi)); wax=max(max(wi));win=min(min(wi));dif=(wax-win)/20;
colorbar('FontSize',15); caxis([win-dif,wax+dif])
shading interp;colormap hot
%title('$$\tilde{\omega}$$','Interpreter','latex')
ylabel('x','FontSize',15)
set(gca,'XTick',0:12:60);xlabel('y','FontSize',15)
set(gca,'FontSize',15)
drawnow
print(gcf,'-dpng',[saldir1,sprintf('%d',0)]);%*dt)]);
    
close

for nt=1:nts
    if (restart==1)
        nt=nt+nts_start;
    else
        nt=nt;
    end

    load([saldir,'dat',sprintf('%4.4d',nt)])
    t=nt*dt;tx=['time=',num2str(t)];
    figure('name',tx)
    colormap jet;

    pcolor(yy,xx,squeeze(wi)); colorbar; shading interp;
    title('$$\omega$$','Interpreter','latex')
    ylabel('x')
    xlabel('y')
%     caxis([0 1])
    drawnow
    
    
%     suptitle(tx);
    drawnow
    print(gcf,'-dpng',[saldir1,sprintf('%d',nt)]);%*dt)]);
    
    close
end